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How localized is an extended quantum system ? 
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We elaborate on a geometric characterization of the electromagnetic properties of matter. A 
fundamental complex quantity, zl, is introduced to study the localization properties of extended 
quantum systems. Zl, which allows us to discriminate between conducting and non-conducting 
thermodynamic phases, has an illuminating physical (and geometric) interpretation. Its phase can 
be related to the expectation value of the position operator (and a Berry phase), while its modulus 
is associated with quantum electric polarization fluctuations (and a quantum metric). We also study 
the scaling behavior of zl in the one- dimensional repulsive Hubbard model. 
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Introduction. The continued cross-fertilization be- 
tween topology and physics has generated simple ideas 
that allow one to clarify many subtle issues in condensed 
matter and quantum field theories. Several problems in 
distinct areas of physics can be phrased in geometrical 
terms resulting in a clearer understanding of its structure 
and a more elegant expression of its solution. A purpose 
of this paper is to show that the concept of localization 
and the field of quantum phase transitions belong to that 
class of problems. Physically, what distinguishes an in- 
sulator from a metal, or superconductor, is the way the 
system responds to external electromagnetic probes, e.g., 
in an insulator a steady current cannot flow at zero tem- 
perature. That macroscopic characterization of matter 
translates microscopically into a study of the localization 
properties of their quantum states, which differ in a fun- 
damental way. How localized is an extended quantum 
system? What are the topological differences between 
the states defining a conductor and an insulator? Can 
one establish a criterion for localization? are some of the 
questions we want to address. The ultimate goal is to es- 
tablish the foundations that will allow one to study and 
predict the electromagnetic properties of generic models 
of interacting quantum particles from general principles. 

The seminal work of Kohn on the nature of the in- 
sulating state is a landmark in the quantum theory of 
matter |jj. Kohn was apparently the first to point out 
that the sensitivity of the ground-state energy, Eq, to 
changes in the boundary conditions is a measure of the 
electric static conductivity of the system, i.e. its Drudc 
weight D c . In fact, in the thermodynamic limit (TL), he 
argued that D c is finite for a non-disordered conductor 
(metal or superconductor), while it vanishes for an insu- 
lator. Physically, this last result is a consequence of the 
localization properties of the many-electron wave func- 
tion, which is exponentially localized in the insulating 
case. Almost thirty years later, Scalapino and co-workers 
H suggested a complementary criteria to Kohn's. They 
proposed to study the limiting behavior of the param- 
agnetic current-current correlation function A xx . The 



appropriate long-wavelength and low-frequency limits of 
A xx should distinguish between an ideal metal, super- 
conductor and insulator. Moreover, they speculated that 
if the system supports a gap one would expect that D c 
equals the superfluid weight. 

Recently, Resta [pLpj and the present authors H have 
introduced a complex quantity zl that allowed one to dis- 
criminate between conducting and non-conducting ther- 
modynamic phases; \zl\ — * or 1 in the large system 
size limit, depending upon the system being conducting 
or insulating, respectively, zl has an enlightening physi- 
cal interpretation: its phase is related to the macroscopic 
electric polarization of matter (Berry phase) while its 
modulus provides information on the polarization fluc- 
tuations (related to a quantum metric nfl). 

It is evident that these different indicators are funda- 
mentally inter-related so, in a sense, they are conceptu- 
ally complementary to each other. Then, why another 
criteria for localization. Apart from helping one to un- 
derstand the nature of metal(superconducting)-insulator 
transitions better, some of the motivations are the fol- 
lowing H: a) physical connection to the modern theory 
of macroscopic polarization |qj|; b) appealing geomet- 
ric content; c) sharper distinction between conductors 
and non-conductors (a two value criteria); d) numerically 
more convenient. 

We have also identified and used geometric concepts 
to characterize the internal topology of stable phases of 
models of interacting quantum particles. Each topolog- 
ical quantum number (related to a Berry phase) labels 
a thermodynamic state, and transitions between them 
determine the boundaries of the quantum phase dia- 
gram [n0[. We have examined phase stability and the 
electromagnetic properties of various strongly correlated 
fermion systems. Here, the study of superconducting and 
Mott phase transitions are among the most fundamen- 
tal issues. In a few cases, exact results have helped to 
elucidate the nature of these transitions, but in general 
one has to rely on numerical calculations of finite sys- 



LA-UR-99-XXXX 



February 1, 2008 



submitted to Physica Status Solidi 



tems for which quantities like D c , or any other corre- 
lation function, vary smoothly at the transition. Con- 
sequently, for example, the boundaries between charge- 
density-wave (CDW) or spin-density-wave (SDW) insula- 
tors and metallic phases are difficult to establish using lo- 
calization indicators, while the boundaries are extremely 
sharp when using topological quantum numbers. 

In the present manuscript we plan to review these fun- 
damental ideas from a general perspective, illustrating 
main concepts in the context of Hubbard-like models. 

General Framework. We will be concerned with generic 
quantum TV-body systems (in a.u. fi, = m = e = c=l)in 
the Schrodinger picture, described by the non-relativistic 
Hamiltonian 



It turns out that 
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where, for simplicity (unless otherwise stated), we will 
assume a one-dimensional (Id) system with particle po- 
sition and momentum labeled by Xi and p il respectively. 
The generalization to higher dimensions is straightfor- 
ward. In condensed matter physics we are interested in 
the behavior of extended systems, i.e., their large- N limit 
(N — > oo and volume L — > oo, such that N/L — hq is 
finite). The boundary conditions on the state $ in L 
then become an essential ingredient. It turns out that 
the generalized Bloch boundary conditions provides the 
system's configuration space M. (of volume L) with a 
non-simply connected structure. Therefore, the state \t 
(G C2(M,N)) becomes a section of a U(l) fiber-bundle 
whose base space is M.. In this way, the topologically 
interesting properties of the bundle will be related to the 
localization properties of the physical system. 

Consider a finite iV-particle system and define the com- 
plex quantity 



zl 
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so that the phase of zl is related to the macro- 
scopic electric polarization B while its modulus pro- 
vides information on the quantum polarization fluctua- 
tions ||,|l^,0]. Kudinov using the fluctuation-dissipation 
theorem showed that, in the TL, a((X 2 ) — (X) 2 ) is finite 
for an insulator while it diverges in a conductor UM . 

The analyticity assumption of zl about a = in the 
TL is equivalent to saying that the system has a gap (see 
below). In general, this is true for finite systems, par- 
ticularly in Id. As we will argue now, (barring possible 
pathologies) zl — > when the system is a conductor. Ba- 
sically, the operator e laX shifts all one-particle wave vec- 
tors by —a. Therefore, in Fermi or Luttinger liquids, the 
Fermi surface is shifted and zl should vanish for L — > oo . 
It is also easy to see that e %aX applied to the BCS wave 
function for the superconducting ground state converts it 
into a current carrying superconducting state orthogonal 
to the former in the TL, which also leads to zl — ► 0. In 
contrast to the Drude weight, these results are also valid 
in the presence of disorder. 

In the case of a system which is insulating due to the 
effect of disorder, with correlations playing a secondary 
role, one expects that the ground state is a single Slater 
determinant with the N particles occupying the N en- 
ergetically more favorable localized wave functions. We 
calculate zl in a more general case, in which the ground 
state is a combination of M Slater determinants 
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where a=%l (n = N/L = n/l, (n, £) are integers with with £* \ C J I 2 = \ and thc single-particle orbital I of the 



n/l an irreducible fraction), X = X) 7 =i^j) an d l^o) is 
the ground state (H |* ) = E |* ), (*o|*o) = !)■ Con- 
sidered as a continuous function of a, zl plays the role 
of a characteristic function generating all moments of the 
operator X, and trivially satisfies zl[Q] = 1, I^lHI < 1, 
zl[(x] — z 1[—a]. Note, however, that the operator X is 
not a genuine operator in the Hilbert space bundle de- 
fined above, although its exponential is a legitimate one. 
Therefore, expectation values of arbitrary powers of X 
have only meaning in terms of z^. Assuming analyticity 
in the neighborhood of a = 0, zl can be written in terms 
of cumulants Cj(X) 



j Slater determinant satisfies 



\1>3i(x)\ < A e-l-^'IA, 
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Then, 
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Because of the bound (|g) , the exponential in the matrix 
element can be expanded in powers of a and only the 
first term survives in the TL L 3> 2ir£A (aX <C 1). Thus, 
orthonormality of the one-particle orbitals leads to 



z L = exp 
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This leads to \zl | — ► 1 in the case of only one Slater de- 
terminant, or in the general case provided that a can be 
chosen in such a way that the exponents in Eq. (||) have 
the same value for all determinants. For translationally 
invariant systems, this is the reason for the introduction 
of the factor £ in the definition a = 2n£/L M. 

For insulating systems in the presence of disorder and 
correlations, we are not able to prove that \zl\ — > 1 with- 
out introducing additional hypotheses. 

On the other hand, under the assumption that |^ (/3)) 
is non-degenerate in the interval < (3 < a — 2tt£/L, for 
any correlated insulator 1^:^ | — > 1 in the large system size 
limit. To understand this statement let's first rewrite Zl 
in the following manner 



2 N 

zl = — (*o|P*o(«)) , P = J>, 

2 — 1 



(9) 



an expression that is obtained after simple algebraic ma- 
nipulations, using the fact that the state |\&o(o0) 1S an 
eigenstate of the a-dependent Hamiltonian 



H = H-aP 



(10) 



with eigenvalue E Q = E - a 2 N/2 (H \^ a (a)) = 
Eol^ (a))) . Assuming that this state remains non- 
degenerate for < a < 2ir£/L, one can iterate the equa- 
tion 

|*o(a)> = |*o> -aG[l- |* (a))<*o|]£|*o(a)> - 

(11) 

with G — (Eq — H) . This is the starting point of 
Rayleigh-Schrodinger perturbation theory to any order 
in a. For an insulating state the result is 



where the Drude weight (or charge stiffness) is defined as 
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According to Kohn's criteria lirn£,_>oo D c — in the 
insulating state, thus, Eq. ( |l3| ) is telling us that 
liniL—too \zl I = 1- It is useful to bear in mind that Eq. 
( |l2| ) is only valid if the ground state never becomes de- 
generate with any other state p|. This is an important 
physical assumption. 

In order to understand the geometric content of the 
localization indicator zl, we need to understand the Rie- 
mannian structure of our Hilbert space bundle. Consider 
a set of normalized states {|4 , o(«))}- Let's assume that 
this manifold of quantum states is generated by the ac- 
tion of the group of Galilean transformations, i.e., 



l*o(a)> 
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where a is a vector and d represents the number of gener- 
ators X (usually the dimension). What is the "distance" 
between two of these quantum states? To answer this 
question, we have to define a (gauge invariant) metric 
tensor in such a manifold. 

The problem of establishing a Riemannian structure 
on an arbitrary differentiable manifold of quantum states 
has been addressed in the eighties |l2| . It turns out that 
a meaningful definition for the metric is 



g^(at) = Re(5 /J *o(o:)|9 y *o(a)) - ^(0)^(0) 



(16) 
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where |\l/j) are eigenstates of H, Py = (\tj|P|\|J,-) and 
Eij — Ej —Ej. Notice that since the parameter a ~ 1/L, 
the 0(a 2 ) does not have the usual meaning because the 
matrix elements themselves depend upon L. However, we 
argue that when L 3> £ (£ is a correlation length), in the 
insulating case, this expansion is well-behaved. There- 
fore, up to third order in a (including the normalization 
of |*o (a)) to the same order) 
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where 7 M (a) = — i(*o( a )l<9^*o( a )) is the Berry con- 
nection p] and d^ — d/da^. It can be shown that g^ v 
is a symmetric and positive definite second rank tensor. 
Moreover, what is of interest to us is the fact that the in- 
finitesimal distance is related to the quantum fluctuations 
of the generator X, i.e., the polarization fluctuations 



g„ v (0) = (X fl X v )-(X^)(X v ) , 



(17) 



and the expectation values are evaluated over |* (0)) = 
|*o)- In a sense, the metric structure on the manifold is 
fixed by the quantum fluctuations which determine the 
modulus of zl in the TL. On the other hand, the an- 
tisymmetric tensor fi M „ = Im(9 M * ( a )|^*o( a )) plays 
the role of a curvature, fl^ is a quantity connected to 
the non-dissipative part of the conductance in adiabatic 
charge transport. 

The Berry phase is an example of a general geometric 
concept which finds realization in several physical prob- 
lems. It is the anholonomy associated with the parallel 
transport of a vector state in a certain parameter space. 
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Anholonomy is a geometric concept related to the fail- 
ure of certain variables to return to their original values 
after others, driving them, are cyclically changed. In con- 
densed matter, the charge Berry phase j c is a measure 
of the macroscopic electric polarization in band or Mott 
insulators while the spin Berry phase 7 S represents the 
difference between the electric polarization per spin up 
and down. j c is determined by the phase of zl , through 
the relation 7 C = a(X). In systems with inversion sym- 
metry 7 C and 7 S can attain only two values: or 7r 
(modulo(27r)). Thus, if two thermodynamic phases differ 
in the topological vector 7 = (j c , j s ) this sharp difference 
allows us to unambiguously identify the transition point 
even in finite systems. This topological "order parame- 
ter" was recently used by us to detect metallic (supercon- 
ductor), insulator and conductor-insulator transitions in 
Id lattice fermion models [ [l0| . 

Application to Lattice models. To illustrate the con- 
cepts developed in previous section, we will study the 
localization properties of the well-known Id Hubbard 
model at half- filling (no = 1) |13, whose ground state 
is always insulating. We concentrate on the analysis 
of the scaling properties of zl as a function of the on- 
site repulsion U/t. For the Drude weight, the scaling 
D c ~ \/Le~ L '^ u > has been established by Stafford and 
Millis H). Using Eq. (|l|), the fact that odd terms 
in a in perturbation theory vanish due to the inversion 
symmetry of the model (hin^^oo zl = ±1), and the ad- 
ditional dependence on L (independent of a = 2n£/L) of 
the sums entering Eq. (|l3|), one finds that for large L, 
\zl\ scales as 1 — £/L, where for large U, £ ~ U~ 2 . 




FIG. 1. Scaling behavior of zl and D c . (a) For the ex- 
ponential(parabolic) scaling: £(4) = 3.93583(4.00134) and 
£(8) = 0.842294(0.797996). (b) £(4) = 3.74949(4.0594), 
£(8) = 1.34926(1.4136) and £(10) = 1.06704(1.1171). We 
quote in parentheses their TL values. 

In Fig. la, we display z L for 4 < L < 24, calculated 
with exact diagonalization (L < 12) and density matrix 
renormalization group (DMRG) (L > 12) for U/t = 4 
and U/t — 8. For both values, a parabola in 1/L fits 
very well the results. Moreover, the values of zl extrap- 
olated to the TL are very close to one (c(4)=l. 00979 and 
c(8)=0. 997598). We also notice that in the scaling regime 
an exponential of the form \zl\ — e~^' L fits the data very 



well (c(4)=1.00993 and c(8)=0.99968), but we do not at 
present have any justification for this functional form. 
For completeness, Fig. lb shows D c for L < 12. 

The ability of z L to detect Mott metal-insulator transi- 
tions as a function of a parameter has been demonstrated 
in the extended Hubbard model for infinite on-site repul- 
sion U/t at quarter-filling || and in the Hubbard model 
with correlated hopping at half- filling p0| . In the latter 
case, due to the pseudospin SU{2) symmetry, the transi- 
tion is best determined by the abrupt jump in the charge 
Berry phase j c . This jump coincides with the change of 
sign of z L for L — > 00 (for any finite system, in general 
zl y^ 0, although \zl\ can be negligible in the metallic 
phase) p|. In the Hubbard model with correlated hop- 
ping, to extrapolate z L to the TL, it is better to take 
U — U C (L) fixed rather than U fixed, where U C (L) is the 
value of U for which z L (U c ) = 0. 

Summary. Over thirty years ago, Kohn [BJ recognized 
that the essential microscopic property that distinguishes 
the insulating state of matter is the disconnectedness of 
its ground-state wave function. In the present paper, we 
have shown that one can quantify this property using 
concepts borrowed from topology. In this way, we have 
seen that the phase of zl (that is a Berry phase) provides 
information on the macroscopic electric polarization of 
matter, while its modulus (connected to a quantum met- 
ric) measures the degree of localization in terms of the 
electric polarization fluctuations, z L displays a qualita- 
tively different behavior for conductors than for insula- 
tors, thereby providing a two- value criterion to distin- 
guish between those states of matter. Moreover, we have 
seen that z L can be computed in a simple and very effi- 
cient way, using standard many-body techniques. In par- 
ticular, we have studied the scaling behavior of \zl\ for 
the Id repulsive Hubbard model in the insulating phase, 
and concluded that it is power law in 1/L. We have also 
briefly mentioned the use of topological quantum num- 
bers (the phase of z L being one of them) to determine 
quantum phase diagrams of interacting particle systems. 

We thank J. Eroles and I. Souza for useful discussions, 
J. Gubernatis for a careful reading of the manuscript, 
and K. Hallberg and CD. Batista for their help with the 
DMRG calculations. Wo also thank T. Souza for poiut- 
ing out Ref. Q to us. G.O. acknowledges support from 
DOE. A.A.A. is partially supported by PICT 03-00121- 
02153 of ANPCyT and PIP 4952/96 of CONICET. 

[1] W. Kohn, Phys. Rev. 133, A171 (1964). 

[2] D.J. Scalapino, S.R. White and S.C. Zhang, Phys. Rev. 

Lett. 68, 2830 (1992). 
[3] R. Resta, Phys. Rev. Lett. 80, 1800 (1998). 
[4] R. Resta and S. Sorella, Phys. Rev. Lett. 82, 370 (1999). 
[5] A. A. Aligia and G. Ortiz, Phys. Rev. Lett. 82, 2560 

(1999). 
[6] N. Marzari and D. Vanderbilt, Phys. Rev. B 56, 12847 

(1997). 
[7] I. Souza, T. Wilkens and R.M. Martin (unpublished). 



LA-UR-99-XXXX February 1, 2008 submitted to Physica Status Solidi 

[8] R. D. King-Smith and D. Vanderbilt, Phys. Rev. B 47, [11] E.K. Kudinov, Sov. Phys. Solid State 33, 1299 (1991). 

1651 (1993). [12] J. P. Provost and G. Vallee, Commun. Math. Phys. 76, 
[9] G. Ortiz and R. M. Martin, Phys. Rev. B 49, 14202 289 (1980). 

(1994). [13] C.A. Stafford and A.J. Millis, Phys. Rev. B 48, 1409 
[10] A. A. Aligia, K. H allberg, CD. Batista and G. Ortiz, (1993). 



cond-mat/9903213 and 9910180 



